Rigorous construction of ground state correlations in graphene: 
renormalization of the velocities and Ward Identities 
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We consider the 2D Hubbard model on the honeycomb lattice, as a model for single layer graphene 
with screened Coulomb interactions; at half filling and weak coupling, we construct its ground state 
correlations by a convergent multiscale expansion, rigorously excluding the presence of magnetic 
or superconducting instabilities or the formation of a mass gap. The Fermi velocity, which can be 
written in terms of a convergent series expansion, remains close to its non-interacting value and 
turns out to be isotropic. On the contrary, the interaction produces an asymmetry between the 
two components of the charge velocity, in contrast with the predictions based on relativistic or 
continuum approximations. 
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The recent experimental realization of a monoatomic 
graphitic film [1], known as graphene, has elicited an 
enormous interest in the study of the properties of two- 
dimensional electron systems on the honeycomb lattice, 
which is the typical underlying structure displayed by 
single-layer graphene sheets. Graphene is quite differ- 
ent from most conventional quasi-two dimensional elec- 
tron gases, because of the peculiar quasi-particles disper- 
sion relation, which closely resembles the one of massless 
Dirac fermions in 2 -f 1 dimensions [2]. Already in the 
absence of interactions, the system displays highly un- 
usual features, such as the anomalous dependence of the 
cyclotron mass on the electronic density, an anomalous 
integer quantum Hall effect, and the insensitivity to lo- 
calization effects generated by disorder. In real systems, 
despite the unavoidable presence of electron-electron in- 
teractions, such consequences of the relativistic-like dis- 
persion relation have been experimentally verified [3, 4]; 
the observation based on angle-resolved photoemission 
spectroscopy [5, 6] are compatible with linear dispersion 
relation and isotropic Fermi velocity. 

From the theoretical point of view, the effect of the in- 
teraction has been analyzed by mean field and Renormal- 
ization Group methods [7-16], but such results are only 
approximate and some non perturbative effects, such as 
a mass generation, have not been excluded so far and 
are still subject of an active debate [17-19]. It would be 
desirable to substantiate the predictions of theoretical 
analyses by rigorous results and exact solutions, which 
would clearly be useful bench marks, but these are rare. 
There are very few rigorous results about the structure 
of the ground state of the Hubbard model in two or more 
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dimensions, among which we would like to mention the 
results in [20] , guaranteeing the uniqueness of the ground 
state and the vanishing of its total spin. However, as far 
as we know, so far no results were proved about existence 
or non existence of long range order and about the long 
distance behavior of correlation in the Hubbard model 
on the honeycomb or other lattices. 

Two theorems reporting the first rigorous construction 
of the ground state correlations in the 2D Hubbard model 
on the honeycomb lattice at half filling and weak cou- 
pling are stated here. Our results exclude the presence 
of magnetic or superconducting instabilities and the for- 
mation of a mass gap. The interaction changes by a finite 
amount the wave function renormalization and the Fermi 
velocity, which remains isotropic, even though the model 
breaks the invariance under 90° degree rotations; on the 
contrary, we predict that the charge velocity develops an 
asymmetry between the two components, an effect that 
is in principle accessible to experiments. 

The analysis is based on the methods of constructive 
renormalization, which have already proven to be quite 
effective in analyzing ID interacting fermionic systems 
in their ground state [22] and 2D systems up to expo- 
nentially small temperatures or in the presence of a non- 
symmetric Fermi surfaces [23-25]. While constructive 
renormalization is based on the ideas of Renormalization 
Group (RG) , the way in which it implements these ideas 
is slightly different with respect to other more standard 
schemes, the main advantages being that the resulting 
method is: (i) exact, in the sense that it does not need 
any relativistic approximation or continuum limit and it 
allow us to keep the full lattice structure of the problem; 
(ii) non-perturbative, in the sense that it involves expan- 
sions whose convergence can be mathematically proved. 

The Hamiltonian of the 2D Hubbard model on the hon- 
eycomb lattice at half filling in second quantized form is 
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given by: 
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where A is a periodic triangular lattice with basis ai,2 = 
|(±\/3, 3) and the nearest neighbor vectors 5i are defined 

as Si = (0,1), 62,3 = ^(±\/3, — !)• The creation and 
annihilation fermionic operators with spin index a =ti) 
ai , 6^ -, , satisfy periodic boundary conditions in x. 

The choice of the interaction is done only for definite- 
ness (it is the simplest one for which the anisotropy of 
the charge velocity is visible at first order in renormal- 
ized perturbation theory) but our results are valid for a 
generic short range density-density interaction. 

We introduce the two component fermionic 
^± ■ ' ' 

^ X,( 



operators ib"^ = f ai , 6^ ? ) and ib^ 
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If 



.^,|AHooTr{e-'3H*T{-}}/Tr{e-^'H'^}, with T the 
fermionic time ordering operator, the zero temper- 
ature 2n-point Schwinger functions arc defined as 
(nr=i V'xi,CTi'0y. ct')- ^'^^ interacting U = case, 

the Fourier transform of the 2-point Schwinger function 
is given by 
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With V-t^ = /_T/2rfa;oEjeAe^*"'''^ 
X]i=i s''^^'^*"'^^'- 'S'o(k) is singular at the Fermi points 
p| = {0,Pf), where ^± = (±^, f ). Close to pp, 

v(k' +pp ) ~ (3/2)(±fc5^ +ifc2), so that the free Schwinger 
function is asymptotically the same as the one of masslcss 
Dirac fermions in 2 + 1 dimension. 

The density operator is defined as pp = (/3|A|)^^ 
X^k (7 V'k o-V'k-p o- 8.nd the definition of the current, for 
?7 = 0, is obtained from the equation rfpx/ dxo = [Ha , Px] : 
in fact, the latter, for small p, assumes the form of a con- 
tinuity equation provided that the current is chosen as 

ip,i = (/3|A|)-iEk,^^k,a<^i^k-p,a' « = 1'2, with £71, (72 
the first two Pauli matrices. The continuity equation im- 
plies the validity of an asymptotic Ward Identity: defin- 
ing k' = k — p^, if k', p, k' — p are small and of the same 
order of magnitude [26], 

((ipoPp ± (3/2)piip,i + (3/2)p2ip,2); V^k.^V^i^-p,.) - 

^ [(v^k,.^^.) - i^'^-pJLpJ] ■ (3) 

When the interaction is present, the Schwinger func- 
tions are not exactly computable anymore; however, 
quite remarkably, they can be computed in terms of a 



convergent rcnormalizcd perturbative series, and their 
long distance asymptotic properties can be rigorously de- 
rived, as summarized by the following theorem. 

THEOREM 1: There exists a constant Uq > Q such 
that, if \U\ < Uo, the specific ground state energy and 
the zero tem,perature Schwinger functions of model (1) 
are analytic functions of U. The Fourier transform of 
the 2-point Schwinger function S(k) is singular only at 
k = and, close to the singularity, it can be written as 
[26], 



5(k)^ 
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where Z = Z(U) and vp = vf{U) are analytic functions 
of U, such that Z = 1 + 0{U^) and vp = 3/2 + bU + 
0{U^), with 



dk 



B, 2|Bi| \v{k)\ 
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(5) 
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and Bi the first Brillouin zone. 

The result summarized in Theorem 1 says that the in- 
teraction does not qualitatively change the asymptotic 
behavior of the 2-point Schwinger function close to the 
Fermi points; the effect of the interaction is essentially to 
change by a finite amount the wave function renormaliza- 
tion and the Fermi velocity. This implies that the inter- 
acting correlations decay as fast as in the non-interacting 
case and, therefore, the presence of quantum instabilities 
in the ground state, such as Neel or superconducting long 
range order, is rigorously excluded at half filling and weak 
coupling, together with the possibility of a mass genera- 
tion. 

Note also that, in the presence of the interaction, the 
Fermi velocity remains the same in the two coordinate 
directions even though the model does not display 90° 
discrete rotational symmetry, but rather a 120° rota- 
tional symmetry; such a remarkable property can be 
easily checked at first order (replacing dp^ by —idp^ in 
Eq.(5), the same result is found); for a proof at all orders 
in the convergent expansion for vp, see below. 

THEOREM 2: For \U\ < Uo, if k', p, k' - p are small 
and of the same order of magnitude [26], then 

{{ipopp ± t^lPljp,! + V2P23p,2) ; i>U,a'^U-p,a) - 

^ - fep,.^^-p,.)] , (6) 

where z = z{U) = 1 + 0{U'^) and the charge velocity 
wi,2 — 3/2 + 0{U) are analytic in U and vi — V2 = aU + 
0{U'^), with 
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Theorem 2 says that, in the presence of interac- 
tions, a new Ward Identity, different from the non- 
interacting one. is verified, tlic main difference with re- 
spect to Eq.(3) being that the charge velocity {vi,V2) is 
interaction-dependent and different from the Fermi ve- 
locity. Remarkably, its two components are different: 
this anisotropy is related to the presence of the lattice, 
that is to the irrelevant terms in a RG sense, which are 
not negleted in our exact scheme. In fact, if we replace 
v{k'+Pp) by its asymptotic expression {3/2){±k[ + ik2), 
the anisotropy coefficient defined by Eq.(7) vanishes ex- 
actly (and so do the higher order corrections). 

We now sketch the proof of the two theorems above 
(for a detailed proof we refer to [21]). The starting 
point is the well-known representation of the ground 
state energy in terms of a Grassman functional integral: 
eo = lim^jAHoo(/3|A|)"^log/P(rfV')exp{V(V')}, where 
P{dijj) is the Grassman gaussian integration with prop- 
agator 5*0 (k), see Eq.(2), and V{tp) is the quartic inter- 
action Eq.(l). One can compute eo by expanding the 
exponential exp{V(V')} in Taylor series in U and naively 
integrating term by term the Grassmann monomials, us- 
ing the Wick rule; however, by such procedure, it is very 
difficult to take into account the cancellations present 
in the perturbative series. The bounds obtained by this 
"simple" procedure are non-uniform in /3, and they do 
not allow one to take the thermodynamic and zero tem- 
perature limits. Therefore, we set up an iterative pro- 
cedure for the computation of eo, based on (Wilsonian) 
Rcnormalization Group (RG) and involving non trivial 
resummations of the perturbative series. 

The first step is to decompose the propagator S'o(k) as 
a sum of propagators supported close to the two Fermi 
points and more and more singular in the infrared re- 
gion, labeled by a quasi particle index uj = ± (labelling 
the Fermi points) and by an integer ft, < 0, so that 

'S'o(k) = YX<i ~ Pf)' "^^^^ 9^^'' supported on 

quasi-momenta of scale 2^ and, on the support, of size 
ll^af^ll ~ 2"''. At this point, we compute eo by iter- 
atively integrating the propagators g^^\g^~^\... After 
each integration step we rewrite 

(8) 

where P{d'>p^'^^) is the Grassmanian quadratic integra- 
tion with propagator given by 

„(<M^^/^ _ Xh(k') ( -iko ChiTk'i + ik'2)\^ 
^± Zh [ch{Tk[-ik'2) -iko ) 

(9) 

where Xh C*') is a smooth compact support function 

nonvanishing only for |k'| < 2'*; V'^'*^ is the effective 
potential, a sum of monomials of arbitrary order, with 
kernels that are analytic functions of U: analyticity is 
a very non trivial property obtained exploiting anticom- 
mutativity properties of Grassmann variables, via Gram 



inequality for determinants. The scaling dimensions of 
the kernels with external lines are equal to 3 — rie, 
modulo an additional dimensional gain, following from 
the fact that all kernels with > 4 external lines are irrel- 
evant in a RG sense (see [21], Theorem 2). The kernels 

W2^^{k') with rie = 2, which are linearly relevant, can be 
inserted step by step in the gaussian integration, thanks 
to the fact that they have the same reality /symmetry 
properties as the free quadratic action: in particular, it 
is found that W^^\q) = and 

^ ^ ^ \5h{Tk[-ik'2) -izhko J 

(10) 

for suitable real constants Zh,Sh- Note that ^k'l and ik2 
are multiplied by the same constant Sh, which is quite 
remarkable; see [21], Lemma 2, for a proof. Iterating 
the procedure above, wc find recursive equations for 
and Cft; in the h ^ —oo limit, the two running cou- 
pling constants converge to values Z-^o = Z = Z{U) 
and c_oo ^ vp = vp{U), which arc close to their un- 
perturbed values and are analytic in U (again, thanks to 
the fact that all sub-diagrams with rig > 4 are irrelevant 
in a RG sense). This completes the discussion concern- 
ing analyticity of eo. A similar discussion allows us to 
prove analyticity of the Schwinger functions and Eq.(4), 
see [21], Sec.III.D. 

Let us now discuss a sketch of the proof of Theorem 
2. We perform a multiscale analysis similar to the one 
sketched above, with V{iIj) replaced by V(V') — B{(j),J), 

withB(0, J) = (<^+,^-)+(V'+,r)+E'=o(^M..?M)- <i>i. 

two external Grassmann fields, ./x./j, an external com- 
muting field, and jx,/t the current (here jx,o = Px, with 
p the density, see the lines preceding Eq.(3)). 

The iterative integration procedure described above, 
in this case, produces, besides the effective potential 
V^^\ new terms involving the external fields. In par- 
ticular, at scale h, the effective source term is given 

by E^,^'4ti'('^/''4''^ ^it^ = '^z'-/''^/'' by the 

discrete invariance symmetries of the model, see [21], 
Lemma 1. A crucial remark is that, while in a rela- 
tivistic QFT zj^'^ is /i-independent, here it is not, the 
relativistic symmetry being broken by the presence of 
the underlying lattice (i.e., by the irrelevant terms in the 
fermionic action). We find that, in the limit h ~oo, 

zl^^ ^ Z^ = Z^{U), which are analytic functions of U, 
with = 1 + 0(C/2) and Z2 - Zx = (2a/3)C/ + 0{U^), 
with a given in Eq.(7). 

Using the fact that all kernels with four or more ex- 
ternal lines are irrelevant, we find that, for k',k' — p, p 
small, 

(ip,/^; V^k-p..^.) - z^si\,)a^s{k - p) , = 0, 1, 2 

(11) 

with (To = 1 and S'(k) the interacting 2-point Schwinger 

function in Eq.(4). The combination S'(k)(T^S'(k -I- p) is 
asymptotically the same as the vertex of a relativistic 



4 



FIG. 1: The two graphs contributing at first order to Z^^, 
/u = 1,2. The full line corresponds to a diagonal propagator 
and the dotted line to an ofl^-diagonal one. Note that the 

two graphs appear with a different sign, depending on the 
value of 11. The second graph is vanishing in the continuum 
approximation v{k) ~ {3/2){±k[ +1^2). 

QFT with wave function rcnormalization Z and "speed 
of light" equal to vf; therefore, it satisfies a relativistic 
WI, relating it to the derivative of 5(k): 

S{k) [ipo + VF{±Picri + P20-2)] ^(k - p) ~ 

^l[5(k)-5(k-p)] . (12) 



Combining (11) and (12) we get Eq.(6) with vi^2 = 
vfZq/Zi^2 and z = Zq/Z. 

In conclusion, we rigorously constructed for the first 
time the ground state, including correlations, of single 
layer graphene with nearest neighbor electron-electron 
interactions; we proved analyticity in U of the theory 
and rigorously excluded quantum instabilities and the 
presence of mass gaps. While in the absence of interac- 
tion the Fermi and current velocity in the Ward Identity 
are the same, they become different in the presence of 
interaction. We proved that the Fermi velocity remains 
isotropic, while the charge velocity develops an asymme- 
try between the two components, a previously unnoticed 
effect of the interaction, which may be detected in future 
experiments. The methods involved are those of con- 
structive rcnormalization, which we believe to be useful 
also for the case of long range (unscreened Coulomb) in- 
teractions, which is currently under investigation. 

We thank G. Benfatto for many useful discussions. 
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